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A systematic set of microwave measurements of the surface impedance [Z s = R s + iX s ) of 
Y BaiGu-sO%.m single crystals (called YBCO/BZO) grown in BaZrOz crucibles reveal new prop- 
erties that are not directly seen in similar measurements on other YBCO samples. The complex 
conductivity a = a\ — ioi obtained as a function of temperature (T) from the surface impedance 
data shows two key features : (1) A new conductivity peak in <Ji(T) around 80K in addition to 
peaks at 30K and 92K, and, (2) extra pairing below 65K in addition to the onset of pairing below 
the bulk T c of 93.47i" as revealed in 02 (T). These features are present in all 3 Y BCO / BZO crystals 
measured and are absent in YBCO crystals grown by other methods. These results show that in 
addition to pairing at T c — 93.4K, an additional pairing channel opens up at ( ~ 6574"). 

High pressure oxygenation of one of the crystals still yields the same results, and shows that 
the data cannot be due to unwanted macroscopic segregation of O -deficient regions. Systematics 
on three single crystals show that the height of the quasiparticle conductivity peak at 80K in the 
superconducting state is correlated with the inelastic scattering rate in the normal state. Close to 
T c , o~2(T) ~ (T c — T), indicating a mean-field behavior and inconsistent with 3DXY fluctuations. 

A single complex order parameter cannot describe these data, and the results suggest that at least 
two superconducting components with corresponding pairing temperature scales (Ta ~ 6574" and 
Tb = T c = 93.4JC" ) are required. Comparison to model calculations considering various decoupled 
two-component scenarios ( A + B — d + s, s + d, d + d) are presented. The comparison shows 
that the experimental data do not distinguish between these various scenarios, however the data do 
require that one of the components be an order parameter (OP) with nodes in the gap, such as a 
d— wave OP. Fit parameters to the calculations using the different scenarios are presented. These 
components are naturally present in all YBCO samples, however impurities appear to suppress the 
pair density of the low temperature A component and lead to greatly enhanced scattering of the high 
temperature B component. Overall, our results strongly suggest the presence of multiple pairing 
temperature scales and energies in Y Ba^Gu^On.^. 



I. INTRODUCTION 



The order parameter of high temperature supercon- 
ductors has been extensively studied recently, and a con- 
sensus seems to have emerged in favor of a d-wave order 
parameter (OP) However there are some notable 
indications (see ref. (^|| for a summary) which suggest 
that, particularly in the most widely studied material, 
YBa.2Cu3O6.g5, a pure d-wave OP may not occur and 
that there are indications of a multi-component OP (eg. 
s + d). 

Material purity is crucial for studies of fundamental 
physical properties not obscured by impurity-related ar- 
tifacts. This fact has been validated time and again in 
experiments on oxide superconductors as well as other 
materials. It is now accepted that improvement in mate- 
rial quality has often resulted in a better understanding 
of the physical properties in solid state systems. 

The recent growth of Y Ba^Cu^Oi-s single crystals 
in BaZrOz(BZO) crucibles has ushered in a new gen- 
eration of ultra- pure samples B. This growth method 



avoids the critical problem of crucible corrosion and leads 
to single crystals with extremely clean surfaces and purity 
exceeding 99.995%. In contrast, the best crystals grown 
in conventional crucibles like Au and yttria-stabilized- 
zirconia (YSZ) have final reported purities of 99.5 — 
99.95% §,§ A number of experiments on YBCO/BZO 
crystals have revealed new features which are either 
greatly suppressed or not present in YBCO/YSZ sam- 
ples. We have recently observed novel features in the 
microwave surface impedance of YBCO/BZO which 
clearly indicates the presence of two superconducting 
components 0. 

In this paper, we present further evidence on the sys- 
tematics in three YBCO/BZO crystals. The results ob- 
tained confirm our earlier findings 0, viz. : (1) A new 
conductivity peak in u\ (T) around 80K which is corre- 
lated with the normal state scattering rate and (2) Extra 
pairing below 65K in addition to onset of pairing be- 
low the bulk T c of 93. 47^ as revealed in 02 (T). A single 
complex order parameter cannot describe these data and 
the data indicate the presence of at least two pairing pro- 
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cesses, leading to at least two OP components. Model 
calculations considering various simple two-component 
scenarios (decoupled d + s and d + d) and are presented. 
We show that the two-component picture can provide a 
quantitative description of the data, whereas a single or- 
der parameter (of any symmetry) does not. In addition 
the two-component picture can provide a description of 
the data on the earlier generation of YBCO/YSZ crys- 
tals also. Analysis of the data close to T c is not consis- 
tent with 3DXY fluctuations and displays a mean-field 
behavior. Systematics on three single crystals show a 
correlation between the inelastic scattering rate in the 
normal state with the quasiparticle conductivity peak at 
80K in the superconducting state. We show that while a 
T-dependent scattering rate is required to quantitatively 
describe the conductivity peaks, the two-component sce- 
nario provides a natural explanation for the location of 
the 30K and 80K peaks in <Ji(T) as arising from pair- 
ing at 65-fT and 93K respectively, in contrast with ear- 
lier interpretation of YBCO/YSZ data which associated 
the 30K peak with pairing at 93K. Overall, our re- 
sults show that a pure d — wave state does not occur 
in Y Ba2Cu30e.g5 and strongly suggest the presence of 
multiple pairing energies in FBa2Cu3O6.95.We note that 
several microscopic pairing scenarios are consistent with 
the conclusions of this paper. 



feasible with the new Y BCO / BZO crystals. 

(2) "Fishtail effect" in magnetization eliminated in 
high pressure oxygen annealed Y BCO7.Q/ BZO and op- 
timally doped Y BCO I BZO samples @. 

(3) Schottky contribution to specific heat suppressed in 
YBCO7.0/ BZO indicating the total absence of magnetic 
moments 

(4) The microwave conductivities o~\ (T) and 02 (T) ex- 
hibit two distinct features not consistent with a single 
superconducting order parameter 

(5) The vortex lattice imaged with STM shows two dif- 
ferent regions which is either representative of two super- 
conducting components or two types of ordered oxygen 
clusters fll5|| . 

(6) Clear observation of a direct first order melting 
transition in YBCOr.o/BZO from a vortex lattice to liq- 
uid without an intervening glassy state fl6[| . 

(7) Evidence of extremely low J c in YBCO 7 . /BZO 
crystals indicative of very low pinning. 

In an attempt to bring forth the essential differences, 
we have presented a comparison of some material and 
physical properties of YBCO/YSZ and Y BCO /BZO 
crystals in a tabular form in Table I. 

II. EXPERIMENTAL RESULTS 



A. Properties of Y BCO /BZO crystals 

One of the causes for the presence of impurities in 
Y Ba-iCuiOi-s single crystals is the random substitution 
at the Cu chain sites by trace amounts of crucible con- 
stituents like Au during the melt growth process. This 
results in local variation of oxygen vacancy distribution, 
introduction of magnetic moments and other local de- 
fects |p|-|lo|. Although the overall T c and sharpness of 
the superconducting transition may not be affected by 
a combination of these elements associated with impuri- 
ties, important features of the superconducting ground 
state like the order parameter symmetry, scattering, su- 
pcrfiuid density etc. are likely to be influenced. It is 
precisely these local impurities which are eliminated in 
(BZO) grown YBa 2 Cu 3 7 - S (YBCO) crystals thus pro- 
viding an opportunity to probe the intrinsic ground state 
properties free from defects. Elimination of the metallic 
impurities also leaves oxygen stoichiometry as the only 
variable which needs to be controlled M. 

It is important to emphasize at this juncture that sev- 
eral new results have been obtained on these new genera- 
tion crystals by a number of experimenters using a variety 
of probes like thermal, magnetic and electrodynamic re- 
sponse, and have led to a clearer picture of the nature of 
superconductivity in YBa2Cu^OQ.^. These results are : 

(1) Clear imaging of the flux lattice of YBCO using 
low temperature STM 1 1 1 . Thus far this has only been 



The microwave measurements were carried out in a 
\{)GHz Nb cavity using a "hot finger" technique Jl7| . 
We measure the surface impedance Z s = R s + iX s and 
penetration depth A(= X s / iiquj) as functions of temper- 
ature T, from which we extract the complex conductiv- 
ity a s = <7i — io~2- This method has been extensively 
validated by a variety of measurements on cuprates and 
borocarbide superconductors [^8| , p^| . 

Three single crystals (labeled AE103, AE105 and 
AE180 and typically 1.3 x 1.3 x 0.1mm 3 in size) from 
different batches grown in BZO crucible and one crys- 
tal (labeled AEXX) of comparable dimensions grown in 
YSZ crucible were measured to study the systematics in 
the different samples. While standard oxygen annealing 
procedures at atmospheric pressure were followed to ob- 
tain optimally doped crystals with oxygen stoichiometry 
around 6 . 95 in AE103, AE105 and AEXX, the AE180 
sample was annealed at a pressure of 100 bars for 20 
hours. Note that this high pressure annealing results in 
elimination of the formation of oxygen vacancy clusters 
and as a consequence, the "fishtail" anomaly in mag- 
netization is suppressed. The crystals (AE103, AE105 
and AE180) grown in BZO had T c = 93AK and AEXX 
had T c = 92AK. All four crystals exhibited very narrow 
transitions in R S (T) at 10 GHz of < 0.3K. 

The temperature dependence of surface resistance 
i? s (T) and change in penetration depth AA(T) were pre- 
sented for AE103 and AEXX crystals and discussed in 
detail in an earlier publication M. In this paper, we 
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show that the data for AE105 and AE180 crystals also 
show the new features seen in the AE103 single crys- 
tal. The data for all the three high purity YBCO/BZO 
crystals clearly reveal a new bump in the vicinity of 60K 
which is not seen in YBCO/YSZ. Also the estimated 
London penetration depth A(0) for YBCO/BZO sam- 
ples are ~ 1000A which is lower than the value of 1400A 
estimated for the YBCO/YSZ crystal. We obtain the 
absolute value of X s by equating R n = X n above T c . 

Since the absolute values of R S (T) and X S (T) are 
now known, we can plot the amplitude I Z S (T) \ = 
^R S (T) 2 + X S {T) 2 as a function of T. Fig.| shows such 
a plot for the three YBCO/BZO samples along with 
YBCO/YSZ for comparison. The virtue of plotting 
Z S (T) | is that one can clearly see the lower A(0) in 
YBCO/BZO which effectively translates to a higher su- 
perfluid density (n s /m) ~ 1/A 2 . 

In the superconducting state, the complex surface 
impedance is related to the complex conductivity and 
penetration depth: Z s (= R s +iX s ) = ^/i^,ouj/(ai — io^). 
Since the absolute values of R s and X s are known, it is 
possible to extract the complex conductivities: 
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These quantities are important because they enable 
comparison with microscopic theories. The pair conduc- 
tivity <72 is a representation of the superfiuid density 
°2 = (/io^A 2 (T)) _1 and is a convenient way of prob- 
ing how the condensate builds up below the supercon- 
ducting transition. The quasiparticle conductivity u\ is 
determined by the quasiparticle density as well as the 
quasiparticle scattering time. 

In Fig.^j, the temperature dependence of <Ji(T) and 
cr 2 (T) are plotted for all the four samples, AE103, AE105, 
AE180 and AEXX. Two striking features emerge clearly 
from a comparison of the data: 

(i) <72 (T) shows two distinct regions of variation above 
and below - 60 - 70K in AE103, AE105 and AE180 
samples in comparison with AEXX. Given the fact that 
A(0) is lower in YBCO/BZO this implies enhanced pah- 
conductivity in these samples. 

(ii) The normal conductivity 0i(T) shows a new peak 
at around 80K{~ 0.9T C ) in YBCO/BZO crystals which 
is absent in YBCO/YSZ. This peak is greatly sup- 
pressed in AE105 as seen in Fig.|[ We show later that 
there is a correlation between this peak and the normal 
state inelastic scattering rate of the samples. It is to be 
noted that the normal conductivity peak at ~ 30 — 35K 
and a very sharp peak near T c are present in all the sam- 
ples. 



Note that these two new features are both absent in all 
previously measured crystals. This is evident from the 
data on the YBCO/YSZ sample shown in Fig.|. We 
emphasize that the YBCO/YSZ data is indeed "canon- 
ical" of the samples prepared by previous growth meth- 
ods. This is illustrated in Fig.|| which compares the 
YBCO/YSZ data with that measured by us on an un- 
twinned LuBa 2 Cu30e.g5 sample Q and by the UBC 
group on YBa.2Cu3O6.g5 j27| ]. Note the excellent consis- 
tency of data on all previous samples with each other in 
Fig|| and the systematic differences with the new data 
in Fig.H 



III. ANALYSIS AND DISCUSSION 



A calculation of the high frequency conductivity re- 
quires proper incorporation of scattering effects along 
with self-consistent calculations of the gap and density 
of states. In addition anisotropy effects should also be 
included for the cuprate superconductors. While this has 
been done pl|-|23|] for unconventional superconductors in- 
cluding those with a d-wave OP, these approaches are not 
easily amenable to comparison with experimental data. 

Instead, in order to compare with the present exper- 
imental data we use a simpler "two-fluid" model of the 
form El: 



<j(uj, T) — 01 - ia 2 
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where /„ and f s represent the fractions of normal and 
superfiuid (with /„ + f s = 1), and t is the relaxation 
time for normal electrons. In this model, the normal elec- 
trons have damping with the usual Drude conductivity 
at high frequencies, and the superconducting electrons 
have inertia but no damping. The quasiparticle and pair 
conductivities can be numerically calculated from: 



/. = (1 - fn) 



dE 



dc 



where E = y'e 2 + A 2 (</»,T), and A(<f>,T) is the gap 
parameter, and < .. >= Jl" dej) indicates an angu- 
lar average over </>. The gap parameters are given by 
A(0,T) = A S (T) and A(0,T) = A d {T) cos(2<£) for 
s-wave and d-wave superconductors respectively. For 
A S (T) and A<j(T) we use appropriate mean-field tem- 
perature dependences. The validity of this two-fluid ap- 
proach, which is used frequently in s-wave superconduc- 
tors following Mattis & Bardeen |^5[ , has been discussed 
for d-wave superconductors by Hirschfeld, et.al. p2[. 
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A. Comparison with models 

We now discuss the comparison of the experimental 
data with detailed models which are implemented assum- 
ing various conditions for the gap parameters. 

1. Single d-wave order parameter: Its failure to describe the 
present results 

There is a general consensus that the superconducting 
order parameter symmetry in cuprates is not a simple s- 
wave type as is the case in conventional superconductors. 
A wealth of experimental results indicate unconventional 
pairing and more specifically point to the presence of 
nodes in the gap. A strong candidate which has been 
found to account for most of the properties in cuprates 
is the theory of d-wave superconductivity. A single com- 
plex d-wave order parameter has been proposed with a 
gap of the form: A(T) = Ad(T)cos2cj) where <j> spec- 
ifies the orientation of the two-dimensional momentum 
of the Cooper pairs. This expression leads to two im- 
portant signatures, namely, the order parameter goes to 
zero in certain k directions and also changes sign as one 
goes around the Fermi surface. Angle-resolved photoe- 
mission Q and phase sensitive SQUID experiments ||] 
have shown strong evidence for this behavior. For exper- 
iments which measure quantities averaged over /c-space 
(like the microwave penetration depth) , the consequence 
of the d-wave gap will lead to features in the low tem- 
perature dependence which should be consistent with the 
presence of a finite density of states within the gap. 

A linear low temperature dependence in YBCO sin- 
gle crystals first observed by UBC group [^7| (and sub- 
sequently confirmed by others including ourselves in 
YBCO @H and also Bi : 2212 @) has been claimed 
as prime evidence for nodes in the gap, which applies to d- 
wave symmetry among other possibilities. Rigorous cal- 
culations of the microwave conductivity in the framework 
of d— wave theories have also been done to explain the 
features observed in prior YBCO/YSZ samples p2|j23f . 

In Fig|| a plot of (A 2 (0)/A 2 (T)) vs. reduced tempera- 
ture (t = T/T c ) for YBCO/YSZ single crystals is shown. 
A comparison to detailed d-wave calculations has been 
carried out by us in an earlier publication fl§| l . Compari- 
son shows that although the low temperature dependence 
is reproduced by calculations (not shown) using a weak 
coupling d— wave order parameter, the agreement is not 
very good at temperatures T > 0.5T C , suggesting that 
strong coupling effects are needed. 

For the YBCO/ BZO single crystals, the low temper- 
ature penetration depth X(T) is indeed linear for all 3 
samples measured up to ~ 20 — 25K with a characteris- 
tic slope dX/dT of about A.5A/K, similar to that of the 
YBCO/YSZ samples. However deviation from linear- 
ity is observed above 25K for the YBCO/ BZO crystals 



due to the onset of the broad QOK bump. This is clearly 
evident from the plots of 02 (T) presented in Fig.|| which 
shows a non-monotonic dependence for AE103, AE180 
and AE105 samples when compared to the AEXX curve. 
While a single order parameter may be reconciled at first 
glance for the YBCO/YSZ data [Q, it is impossible to 
do so for the results on the BaZrOz grown crystals 0. 
Instead a two-component model has to be considered as 
is done below. 

B. Comparison to decoupled two-component order 
parameters 

The simplest way in which one can analyze a two- 
component system is to consider two parallel supercon- 
ducting channels and to add their complex conductivities. 
The total conductivity can then be written as: 

o = (J\ - i(7 2 = (&IA + °ib) - i(cr 2 A + &2b) 

For ease of calculation we consider that the two com- 
ponents are decoupled with distinct transition tempera- 
tures, T c a and T c b (> T c a) ■ (We emphasize that in all 
likelihood the components are coupled, and furthermore 
there is only one phase transition at T c b, with T c a more 
of a crossover temperature, as discussed later). Calcu- 
lations of the conductivities can be performed using the 
Mattis-Bardeen formalism introduced earlier. In an ear- 
lier paper, we had presented the qualitative model where 
for ease of calculation we used s— wave order parameters 
for both the A and B components 0. In reality, it is es- 
sential to choose at least one component to be d— wave to 
produce the linear low temperature dependence. Using a 
combination of s and d gap symmetries, it is possible to 
obtain excellent quantitative fits to the observed experi- 
mental data for all the YBCO/ BZO single crystals. 

Fig.^ shows the experimental data for 02 and a± along 
with the fits obtained with the two component model as- 
suming an (s+d) case for the A and B components. Good 
fits can also be obtained for (d+d) and (d+s) cases for 
the two components. The data for the AE103 sample 
only is presented in the figures for clarity. It is possible 
to obtain similar fits for the other YBCO/BZO crystals 
and the fit parameters for all the samples are tabulated in 
Table II. A good agreement between the data and model 
is obvious. Best fits are obtained for the following generic 
choices: 

1. A low temperature component A with T c a ~ 60A" 
and a high temperature component B with T c b ~ 
93K. 

2. A weak coupling d— wave gap with = 2.16fcsT c 
for the the B— component and an s— wave gap with 
Aa = 1.76fcsT c or a weak-coupling d— wave gap 
for the A— component. 
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3. A lower fraction of s— wave (~ 0.3) than the 
d— wave part (~ 0.7). 

4. Temperature dependent scattering rates t^ b = 



r A< B = r AM0) (e a ^-e 

for for <j\. This exponential variation agrees with 
the suggestion that the scattering time increases 
rapidly below the transition temperature. The vari- 
ation of ta and tb with temperature is plotted in 
bottom of Fig.^. The detailed functional form is 
not as important as the key feature of rapid varia- 
tion of ta below 65AT and tb below 93A. 

The two-component model also provides an explana- 
tion of the YBCO/YSZ data as can be seen from FigJ|. 
The fit parameters (Table I) suggest that impurities sup- 
press the order parameter for the A component (a smaller 
ratio 10% of the supcrfluid density n s A as well as a weaker 
gap A J 4(0)/fcsT Cj 4 ~ 1.0 are needed), as well as greatly 
enhancing the scattering rate t b x about 10 times larger 
for the B component. The rapid variation of t a below 
50K is still needed to obtain the 30 A peak in o~i{T) for 
the YBCO/YSZ crystals. 

The occurence of the low temperature normal conduc- 
tivity peak in o~\ at around 30A, even though the su- 
perconducting transition temperatures are in excess of 
90A in YBCO crystals, has long been a puzzling feature. 
An exponential reduction of the quasiparticle scattering 
rate below T c was invoked to account for this. How- 
ever, the present data provides an explanation for the 
location of this peak, since in the present experimen- 
tal results, it is natural to associate the two o~\ peaks 
with the corresponding features in the pair condensate 
density as inferred from a%- Two types of pairing with 
different characteristic energy scales are clearly observed 
in these high quality YBCO/BZO single crystals. In 
the YBCO/YSZ samples, these energy scales are also 
present but their signatures are obscured by impurity 
scattering. 

One aspect of the comparison should be noted. Be- 
cause we have used a model in which the components 
are decoupled, the calculations necessarily show a sharp 
break at around T c a ~ 60A, whereas the data display 
a smooth crossover. This is an indication that the de- 
coupled model is too simplistic, and a coupled model 
is necessary. Indeed such calculations performed ]66| ] 
within a Ginzburg-Landau framework lead to a smooth 
crossover in A _2 (T) in closer agreement with the exper- 
iments. However a full calculation of the conductivities 
requires a microscopic model which still needs to be im- 
plemented. (See also the remarks in the Summary section 
of this paper). 

We note that there have been other discussions of mul- 
tiple components (gaps or quasiparticles) in the cuprate 
superconductors. A two-gap model was proposed earlier 
by Klein et al. pi to describe the data on YBCO thin 



r * A B were used 



films with low cation disorder. In their case, however, low 
and high gaps both with s— wave symmetry and the same 
T c were chosen to calculate the pair conductivity and the 
issue of o~\ was not addressed in detail. Two-component 
behavior in 124 and underdoped 123 compounds have 
also been reported in Raman p0[ and femtosecond opti- 



cal response 31 1. 

Anisotropic penetration depth measurements along the 
a and b— axes reported on untwinned YBCO/YSZ sin- 
gle crystals [^2j. Lower A(0) in the 6-direction suggested 
superconductivity in the Cu — O chains with the onset 
T c same as that of the planes viz. 93A". The similar- 
ity between 1/A„ and 1/A 2 was taken as evidence that 
chain order parameter also has nodes in the gap. How- 
ever, in our present data, additional superfluid response 
clearly turns on below 60 A and if this is due to chains, 
then this would imply that chain superconductivity oc- 
curs at a lower temperature than the planes and would 
also directly contradict the conclusions reached in ref. 
J32|. Recent experiments on anisotropic thermal con- 
ductivity in detwinned YBCO/YSZ single crystals have 
also revealed enhanced superfluid density below ~ 55AT 
which has been interpreted as due to chains ]33| . This 
will be consistent with our observations of the microwave 
response in YBCO/BZO. 



C. Near-Tc behavior: 

We turn our attention now to the transition region 
close to T c . The temperature dependence of the penetra- 
tion depth A(T) is expected to reveal the characteristic 
nature of the superconducting gap as it opens up and also 
possible contributions due to fluctuations. In contrast 
to conventional superconductors, the low dimensional- 
ity and small coherence length in cuprates make them 
good candidates for studying the effect of superconduct- 
ing fluctuations and in determining the universality class 
they belong to. While Ginzburg-Landau theory gener- 
ally describes the region near the transition very well in 
conventional superconductors, it has been argued that in 
YBCO assuming the simplest case of a single complex 
order parameter, it is possible to observe critical behav- 
ior which would give rise to fluctuations corresponding to 
the universality class of the three dimensional (3D) XY 
model. Transport and thermodynamic measurements in 
the presence of a magnetic field have been interpreted in 



terms of this model |3Jj3J 
disputed by Roulin et al. 



However, this has also been 
[ 3j| who argue that fits of the 
specific heat near T c in YBCO single crystals for high 
magnetic fields do not weigh in favor of 3D XY scaling. 

Penetration depth data close to the transition region 
will be influenced by the presence or absence of criti- 
cal fluctuations. Specifically, a plot of [A(0)/A(T)] n as 
a function of T just below T c is a useful indicator 
of the validity of standard mean-field or the 3D XY 
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models. While a mean-field behavior can be deduced 
if [A(0)/A(T)] 2 is linear in T, 3D XY would require 
[A(0)/A(T)] 3 to be linear in T. The latter behavior has 
been observed in penetration depth measurements by the 
UBC group |37| and has been interpreted as consistent 
with the critical behavior of the 3D XY model. It is 
important to note the the choice of the range of temper- 
ature over which one looks for effects due to fluctuations 
is crucial. The dynamical fluctuations in the frequency 
dependent microwave conductivity above T c has been ex- 
amined by Booth et al. p8[ . They arrive at the conclu- 
sion that while the Gaussian fluctuations can account 
for the behavior over a larger range in temperature, just 
1 — 2K above T c , the behavior is dominated by critical 
fluctuations. It should be noted that the temperature 
dependence of [A(0)/A(T)] 2 from kinetic inductance mea- 
surements has shown mean- field dependence |3^| . 

In a bid to investigate the transition region in our 
data on YBCO/BZO and YBCO/YSZ single crystals, 
we have plotted the penetration depth for all four crys- 
tals as a log-log plot against the reduced temperature 
[(T c - T)/T c ] shown in the top panel of Fig.g. It is clear 
that all the data show a linear behavior (with a nega- 
tive slope) with typical values of the slope lying around 
0.5. The appropriate A(0) and T c values estimated from 
our R s and X s data for the three YBCO/BZO samples 
(AE103, AE180, AE105) and the YBCO/YSZ sample 
(AEXX) were taken to obtain the individual curves. Ex- 
cept the data for AE180 which has a slightly higher slope 
(~ 0.6), the other data form a set of parallel lines. This 
behavior is consistent with mean-field variation of the or- 
der parameter below T c . To illustrate the point further, 
a plot of [A(0)/A(T)] 2 vs T is presented in the bottom 
panel of Fig.^3] for one YBCO/BZO crystal. The linear 
variation is quite obvious and is clearly different from the 
behavior reported for YBCO single crystals in ref. [p7| . 
The inset shows the same data plotted as [A(0)/A(T)] 3 
vs T which is non-linear and rules out the possibility 
of 3D XY. The smooth variation of the [A(0)/A(T)] 3 
vs T curve also indicates against the possibility of in- 
terpretation as a crossover behavior from 3D XY very 
near T c to a mean-field like variation at lower tempera- 
tures. Our conclusion in this context would be that the 
near-T c behavior in all YBCO crystals is governed by 
standard mean-field expression. It is interesting to note 
that this universal behavior may also be indicative of a 
scenario against the possibility of a single order param- 
eter. Absence of 3D XY correlations is not surprising if 
one considers unconventional pairing states with mixed 
symmetry. 

D. Normal state scattering rate 

Assuming a local j — E relation in the skin depth 
limit, the normal state surface resistance can be writ- 



ten as R n — yJu>HQp n /2 where the microwave normal 
state resistivity is expected to be the same as the DC 
resistivity and p n — 2T p \(0) 2 . Here A(0) is the London 
penetration depth and T is the normal state scattering 
rate. Taking p n = Po + jT the resistivity values can 
be translated into the inelastic scattering rates given by 
r = 7 T/2 Mo A 2 (0). 

In Fig.^, we have plotted the normal state inelas- 
tic scattering rate as a function of temperature for 
lOOif < T < 200K for all the YBCO/BZO samples 
along with the data for YBCO/YSZ crystal. Among 
the YBCO/BZO single crystals, the T values as well as 
the slopes for AE103 and AE180 over the entire temper- 
ature range are lower than that for the AE105 sample. 
The YBCO/YSZ AEXX sample has the largest T in 
this batch of crystals whose microwave surface impedance 
were measured in an identical manner. A comparison 
of this data with the <j\ (T) data of Fig.|] shows a clear 
correlation between the normal state scattering rate and 
the new 80K quasiparticle conductivity peak. This peak 
is prominent only in AE180 and AE103 samples which 
have a lower scattering rate, is suppressed in AE105 
which shows a larger T than the other two YBCO/BZO 
samples and practically absent in AEXX which has the 
largest scattering rate. This remarkable correlation indi- 
cates the sensitive nature of the new conductivity peak 
to very small changes in impurity concentration and/or 
oxygen stoichiometry. Further studies in YBCO/BZO 
crystals with transition-metal doping and variation of 
oxygen stoichiometry are needed to explore the correla- 
tion between the electronic properties in the normal state 
and the quasiparticle contribution in the superconduct- 
ing state which is suggested by our present data. 

IV. SUMMARY OF MICROWAVE PROPERTIES 
OF YBCO/BZO AND YBCO/YSZ CRYSTALS 

We summarize below some of the microwave properties 
exhibited by the YBCO/BZO crystals along with com- 
parison of results on the previous class of YBCO/YSZ 
samples. The key observations are: 

1. The linear behavior of the low temperature pene- 
tration depth A(T) cx T appears to be very robust 
and is observed in all samples with essentially the 
same slope. A(T) rises linearly with T upto 0.25T C 
in YBCO/BZO single crystals with a character- 
istic slope (dX/dT) of about 5A/K. This linear 
variation is in agreement with the observations on 
YBCO/YSZ crystals too and has been ascribed to 
the presence of nodes in the gap, and overall is con- 
sidered a strong evidence for d— wave order parame- 
ter [|7). However, unlike the data in YBCO/YSZ, 
the full temperature dependence for YBCO/BZO 
display a non-monotonic behavior which points to 
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the presence of an extra component in addition to 
a d— wave part. 

A linear behavior of A(T) is also seen in good qual- 
ity BSCCO single crystals A T 2 behav- 
ior which has also been reported in YBCO thin 
films |6^| , has been considered a consequence of in- 
homogeneties and additional scattering caused by 
impurities and disorder. 

2. The surface resistance R S (T) of YBCO/BZO 
shows a bump at around 2>QK as well as a broad 
bump at higher temperatures. In terms of conduc- 
tivity, this results in two peaks occuring at ~ ZQK 
and ~ 80if . The ratio o\/a n of the peaks are much 
higher and cannot be explained due to coherence 
effects predicted by BCS theory. We have associ- 
ated these two peaks with two different quasiparti- 
cle systems corresponding to components A and B 
as discussed earlier. 

In YBCO/YSZ, only the low temperature peak at 
~ 30K is seen. This was first observed at sub-THz 
f33| and at microwave frequencies |3jj] . To account 
for this feature, Bonn et al. J64|] proposed an rapid 
increase in the inelastic scattering time (r) below 
T c followed by a reduction in carrier density at low 
temperatures. 

3. The interpretation of the conductivity peak as aris- 
ing from an interplay of increasing r and decreas- 
ing n qp as T is lowered has also been used here to 
account for the conductivity peak |33|,|4| . Such an 
interpretation necessarily requires a rapid variation 
of tb just below 93K and of ta below 65K as is ev- 
ident from Fig. I in YBCO/BZO. 

It is to be noted that in BSCCO, our results |ll| 
do not show any conductivity peak in o~\ below T c 
but rather a\/a n rises monotonically as T is low- 
ered. However another study |^5(| claims to see a 
small reduction in <j\ /o~ n at very low temperatures 
suggestive of a conductivity peak. 

4. It is important to note that the two temperature 
scales of93K and 65K are present in YBCO/YSZ 
also. This can be seen from the fact that two com- 
ponent analysis also describes the YBCO/YSZ 
data, as shown in Fig.||, provided a suppressed 
n s A and tb are used. This suggests that in the 
YBCO/YSZ crystals, impurities lead to suppres- 
sion of the pairing density of the A component, as 
well as greatly enhanced scattering of the B com- 
ponent. 

5. A very sharp peak in o\ just below T c at around 
91 - %2K is seen in all (YBCO/BZO as well as 
YBCO/YSZ ) samples, so that there are a to- 
tal of 3 conductivity peaks. This peak near T c 



has been ascribed to fluctuations |40| or inhomo- 
geneities fS5| . In the YBCO/BZO samples this 
peak is extremely sharp and further attests to the 
high sample quality. 

6. The temperature scale of 65K naturally raises 
the possibility of a multi-phase sample with re- 
gions of YBCOe.5- We have considered this issue 
very carefully and several experiments show that 
this is conclusively ruled out. In-situ resisitivity 
meaasurements at high temperatures have shown 
that the crystals can be reversibly oxygenated- 
deoxygenated, and that the diffusion constants 
are well characterized, indicating nothing unusual 
about the oxygenation of these samples ||. Note 
also that the YBCO/YSZ sample AEXX was 
oxygenated using the same procedures but did not 
display the new features in the conductivity. Fur- 
thermore, the similarity of the AE180 data with 
the other samples conclusively shows that even lo- 
cal O vacancies are not responsible for the present 
results. This is because in AE180, high pressure 
oxygenation breaks up O deficient clusters, as evi- 
dence by the elimination of the fishtail anomaly in 
samples prepared by this method. 

7. Our measurements on both YBCO/BZO and 
YBCO/YSZ single crystals are consistent with 
mean-field behavior rather than 3DXY close to 
T c . Near-T c variation of 1/A 2 has been analysed 
in terms of 3DXY scaling by Kamal et al. |37[] in 
their YBCO/YSZ crystals. While a more rigor- 
ous study of the fluctuation contributions below 
and above T c by Booth et al. (3^] indicated that 
3DXY scaling is seen only over a limited temper- 
ature range, kinetic inductance measurements by 
Lemberger et al. iffflj on thin films suggest mean- 
field behavior, in agreement with the conclusions of 
this paper. 

8. It should be noted that the pairing onset below QbK 
in the c^iT) is quite broad, and indicates that this 
is not likely to be a critical point. Thus we believe 
that there is only one transition at T c = 93AK, 
with a crossover at 65K where the additional pair- 
ing channel opens up. Furthermore, although we 
have used a decoupled scenario to analyze the data, 
it is very likely that the OP's are coupled. This 
is confirmed by calculations of the supcrfluid den- 
sity using a simple Ginzburg-Landau free energy 
for two coupled OP, which reproduce the main fea- 
tures of the A~ 2 (T) data p(|. At high tempera- 
tures T c > T > T c a, the symmetry is predomi- 
nantly of type B, although a small A— component 
is also induced due to the coupling. At lower tem- 
peratures T < T c a, the A— component grows so 
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that both components (and hence a mixed symme- 
try) are present. Microscopic calcuations of ther- 
modynamic properties of superconducting lattice 
fermions which also have multiple pairing channels 
also show that the thermodynamic signatures, such 
as specific heat cusps, are weaker in the subdomi- 
nant channels which appear at lower temperatures 
than in the main channel p4| . 

While we have presented a simple phenomenologi- 
cal two-component model in this paper to account 
for the two conductivity peaks in YBCO/BZO, 
microscopic calculations (similar to those based 
on a single d— wave order parameter |2l] - |23| ) 
need to be done for the surface impedance of 
a two-component superconductor. Furthermore, 
the analysis used here assumes phase coherence 
throughout the sample and does not include con- 
tributions arising from phase fluctuations |67|j68[| . 



in determining the possibility of multiple pairing ener- 
gies. The presence of pockets with d and g— wave sym- 
metries has been recently discussed Q. 

(4) Surface states and Time Reversal Symmetry break- 
ing: The possibility of mixed order parameter symmetry 
at the surface leading to breaking time reversal symme- 
try has been suggested |55|,^6|]. Recent observation of 
Andreev bound states at the surface of YBCO p7 58 
strongly suggests this possibility. 

(5) Interlayer Tunneling: In a recent paper, Xiang and 
Wheatley []59| presented a model based on proximity ef- 
fect incorporating a microscopic pair tunneling process 
which couples the Cu — O chains and planes in YBCO , 
to account for the experimental results on YBCO/YSZ 
data from the UBC group Our present results on 
the YBCO/BZO crystals suggests that within the inter- 
layer tunneling model, single electron tunneling processes 
could account for the features observed in our data. 



V. RELATION TO MICROSCOPIC THEORIES 



VI. CONCLUSION 



Models based on any type of mixed order parameter 
symmetry would yield results consistent with our experi- 
mental observations in the YBCO/BZO system. Mixed 
order parameter symmetry and multi-component behav- 
ior of the order parameter has been addressed in several 
theoretical papers which can be broadly classified along 
the following categories: 

(1) s + d,s + id : It has been argued that the or- 
thorhombic structure of YBCO naturally leads to a mix- 
ing of s and d order parameter components J3j. Several 
theories [| 41 - 4^] advocate to this general concept of a 
mixed s + d OP and have successfully developed models 
which can describe the experimental results in YBCO. 
This appears to be consistent with photoemission data 
on overdoped BSCCO also [E6LETJ . Phenomenological 
jyj and microscopic theories |47f| have been considered. 
The work of ref. JilJ indicates that the system likely goes 
from a d— wave state at high temperatures to an s + id 
state at low T. 

(2) Chain-Plane coupling: The presence of Cu — O 
chains in addition to planes raises the possibility of dif- 
ferent superconducting condensates with different pair- 
ing energies residing on chains and planes. This can ef- 
fectively lead to a two-component system exhibiting two 
different gaps Q. Coupling between chain and plane 
bands has been proposed j49|] to account for both the s 
and d characteristics displayed by YBCO. 

(3) Multi band models and Fermi pockets: Theoreti- 
cal treatments based on multi band models have been 
presented which considers the interband interactions 
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In summary, we have shown that high quality sin- 
gle crystals of YBCO grown in BaZrO^ crucibles ex- 
hibit features in their microwave properties, which are 
inconsistent with a pure d-wave OP and instead point 
to the occurence of multi-component superconductiv- 
ity. The measurements reveal the presence of two pair- 
ing temperatures corresponding to two superconducting 
(pair/quasiparticle) components in the optimally doped 
compound. This should be borne in mind and careful 
studies should be performed on high quality materials to 
further explore these issues. 
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YBCO/YSZ 


YBCO/BZO 


Crucible: Yttria-stabilized Zirconia (YSZ) 


Crucible: BaZrOz 


Crystal Purity: 99.5 to 99.95 % [5 


Crystal Purity: 99.995 % |8| 


T c = 93.2 to 93.5K 


T c = 93.2 to 93.5^ 


Estimated A(0) ~ 1600 - 2000A 


Estimated A(0) ~ 1000 - 1400A 


Flux lattice not imaged with STM 


Flux lattice imaged with cryogenic STM 11] 


Schottky anomaly present in specific heat 


Schottky anomaly suppressed in high pressure 
oxygenated O6.95 and eliminated in O7 Jli)] 


Large "Fishtail effect" in Magnetization 


"Fishtail" greatly suppressed in high pressure 
oxygenated Oe.95 and eliminated in O7 Jl2| 


Substantial Pinning due to impurities 


Very low pinning evidenced by extremely low critical 
currents and observation of flux lattice melting jl6| 



TABLE I. A comparison of some material and physical properties of Y BCO / BZO and YBCO/Y SZ single crystals. 
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Samples 


AE103 (d + d) 


AE105 (d + d) 


AE180 (d + d) 


Components 


A(0.6), B(0.4) 


A(0.6), B(0.4) 


A(0.48), B(0.52) 


T C (A,B) 


63K, 93K 


59K, 93K 


58K, 93K 


F(a,b)(o) meV 


11, 8 


10, 13 


7, 13 


a 


12, 25 


10, 8 


9, 14 


r (A,s) meV 


0.08, 0.6 


0.08, 1.1 


0.06, 1.2 



Samples 


AE103 (s + d) 


AE105 (s + d) 


AE180 (s + d) 


Components 


A(0.6), B(0.4) 


A(0.6), B(0.4) 


A(0.45), B(0.55) 


T C (A,B) 


63K, 93K 


59K, 93K 


58K, 93K 


F(A,B)(a) meV 


12, 8 


18, 14 


10, 16 


a 


11, 25 


11, 6 


10, 16 


r (A,s) me V 


0.05, 0.8 


0.04, 0.4 


0.02, 1.4 



Samples 


AE103 (d + s) 


AE105 {d + s) 


AE180 (d + s) 


Components 


A(0.78), B(0.22) 


A(0.78), B(0.22) 


A(0.68), B(0.32) 


T C (A,B) 


63K, 93K 


59K, 93K 


58K, 93K 


r<A,s)(o) meV 


20, 7 


12, 9 


10, 10 


a 


11, 21 


11, 10 


9, 14 


r (A,s) me V 


0.24, 0.52 


0.06, 0.01 


0.1, 0.4 



TABLE II. Fit parameters for the three Y BCO / BZO crystals using the several two-component scenarios. For the 
YBCO/YSZ case, fits shown in Fig. fusing the two-component model are generated using the following parameters: Compo- 
nents: A(0.1) and B(0.9); T C A, B = 60, 9274"; F (AiS)(0) = 10, 14; a = 7, 12; V* (A B) = 0.1, 0.01. While the numbers are identical 
for (s + d) and (d + d) cases, (d + s) does not fit the data. 
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FIG. 1. Magnitude of surface impedance amplitude 
| Z 3 (T) | of YBCO/BZO and YBCO/YSZ single crys- 
tals. Filled symbols are used for BZO grown crystals: AE103 
(filled squares), AE180 (filled triangles) and AE105 (filled di- 
amonds). The data for YSZ grown AEXX (open circles) is 
also shown. The plot emphasizes the fact that A(0) for all the 
YBCO/BZO crystals is lower than that of YBCO/YSZ. 
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FIG. 2. Complex conductivities 02 (top panel) and o\ (bot- 
tom panel). The Y BCO / BZO data are shown in filled sym- 
bols and the Y BCO /Y SZ data in open circles. Note the 
extra pairing at around 60K in 02 and the new peak at 80K 
in tri in YBCO/BZO, both absent in YBCO/YSZ data. 
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FIG. 3. Superfluid density of the YBCO/YSZ single crys- 
tal (AEXX) plotted against reduced temperature. Data for 
an untwinned YSZ grown LuBCO crystal and the UBC data 
are also shown for comparison. 
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FIG. 4. <r 2 and tj x data for YBCO/BZO single crystal 
(AE103) along with the fits generated using a decoupled 
two-component model (s + d) [solid line] and (d + d) (dashed 
line), discussed in the text. The bottom panel shows the vari- 
ation of scattering times ta and tb below the T c for A and 
B components for the (s+d) case. Essentially similar fits can 
be obtained for (d + s) case (not shown) . 
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FIG. 5. tT2 and ai data for YBCO/YSZ single crystal 
(AEXX) along with the fits generated using a decoupled 
two-component model (s + d) discussed in the text. The 
bottom panel shows the variation of scattering times ta and 
tb below the T c for A and B components for the (s+d) case. 
The fit is identical with similar parameters for the (d + d) 
case. The (d + s) case does not fit the data. 



FIG. 6. Top panel: Penetration depth near T c shown in 
a log-log scale. The data sets indicate a similar variation in 
all four samples, with a slope of 0.5 indicating mean field 
behavior. Bottom panel: Plot of the superfluid density, 
[A(0)/A(T)] 2 vs T just below T c for AE103 crystal. The solid 
line is a guide to the eye indicating the linear variation con- 
sistent with mean-field behavior. Inset shows [A(0)/A(T)] 3 vs 
T for the same data in which the curvature is clearly evident 
indicating disagreement with the 3D XY model. 
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FIG. 7. The normal state inelastic scattering rate for all 
four samples. The data for AE103 and AE180 have lower 
values and smaller slope in comparison to AEXX and AE105 
samples. 
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